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The problem of linear stability of confined V-flames with arbitrary gas expansion is addressed. 

\ Using the on-shell description of flame dynamics, a general equation governing propagation of dis- 

. turbances of an anchored flame is obtained. This equation is solved analytically for V-flames in 

' high-velocity channel streams. It is demonstrated that dynamics of flame disturbances in this case 

' is controlled by the memory effects associated with vorticity generated by the curved front. The 

_ , perturbation growth rate spectrum is determined, and explicit analytic expressions for the eigen- 

' functions are given. It is found that the piecewise linear V-structure is unstable for all values of the 

' gas expansion coefficient. 

lO ■ PACS numbers: 47.20.-k, 47.32.-y, 82.33.Vx 

(N ; 

I. INTRODUCTION 

' O ■ Among the various types of premixed flame propagation problems, anchored flames hold a special place. On the one 
^ , hand, such flames are relatively easy to realize experimentally; on the other, they look simple enough for theoretical 

■ investigation, because they admit several important simplifications. For instance, open flames anchored by means 
^ [ of a thin rod are often observed to have rectilinear wings (unconfined V-flames). Homogeneity of the upstream 
^ ■ flow, adopted usually as the natural approximation compatible with this piecewise linear flame-front structure, often 

• ^ ' conveys the impression that the problem is easily solvable analytically. It thus represents an excellent laboratory for 
testing our understanding of premixed flame dynamics. 

Despite these promising circumstances there is an apparent lack of theoretical results on V-flame properties. The 
I reason is that these flames are not as simple as they seem. A closer inspection of the flow structure of the idealized 
V-configuration reveals that this pattern is singular: the pressure field turns out to diverge logarithmically near the 
I \ tip of the fiame- front (and also at infinity along the front, in the case of unconfined V-fiames). This is a sign of 
^ ■ incompleteness of the idealized picture, which means that the system anchoring the flame must be explicitly included 
^sD \ iiito consideration. This essential complication necessitates the introduction of a specific inner scale in the problem 
i (in addition to the cutoff wavelength and the channel width) , thereby raising the question as to the influence of this 
new scale on the whole basic pattern. The initial problem is thus naturally divided into two parts: 1) modeling of 
, the anchoring system; this primarily is a stationary analysis, aimed at inferring properties of the system needed to 
' generate a presumed flame pattern, and 2) investigation of the flame dynamics, which first and foremost is a stability 
] analysis of the anchored flame; an important issue in this analysis is its model-dependence, i.e., the extent to which 
its results depend on particularities of the anchoring system. 

The purpose of the present paper is to carry out an analysis of the above-mentioned issues, in the case of a 
confined V-flamc anchored in a high-velocity gas stream. It will be shown that the problem admits a full theoretical 
investigation in this important particular case, and that its results are model-independent in the above sense. It should 
^ j , be mentioned that in contrast to unconfined anchored flames, flames anchored in channels do not exhibit an acute 

■ linear structure, although the piecewise linear front with a uniform upstream flow is still a solution of the governing 
equations. Experiments show that deviations from linearity occur not only in the small regions near the anchor and 
the channel walls, but all along the front P, H, Q. This suggests that the simplest configuration is possibly unstable 
in the confined case. The results of our work fully confirm this conjecture. 

In our investigation, we use the on-shell description of flames developed in Refs. [1, [1, d, 0|- The integro-differential 
equations derived therein provide a non-perturbative description of spontaneous flame dynamics in the most general 
form, i.e., they apply to flames with arbitrary gas expansion and arbitrary jump conditions across the flame front. 
The main advantage of using these equations is that they are closed, in the sense that they involve only quantities 
defined at the flame front. This allows one to avoid explicit solving of the flow equations in the bulk, which is the 
stumbling block of conventional analysis. This approach will be shown to be extendable to the case of anchored flames 
in a simple and natural way. 

The paper is organized as follows. Section|lT]serves to set up the general framework of the on-shell flame description. 
In Sec. Ill Al we formulate the problem and recall the main results of Refs. 0, H, H, 0|- Extension of these results 
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to the case of anchored flames is described in Sec. Ill Bl An analysis of the anchoring system impact on the flame 
structure, carried out in Sec. Ill C[ is used in Sec. HID II to identify boundary conditions for the Hnearized equation 
describing the propagation of disturbances. This equation is derived, in a form suitable for the subsequent analysis, 
in Sec. HID! then solved in Sec. IIIII An important step here is the evaluation of rotational contribution, presented 
in Sec. IIII Al The resulting equation is analyzed in the high-velocity limit in Sec. IIII Bl which allows considerable 
simplifications. In particular, an asymptotic expansion of the main integral operator J{ is constructed in Sec. IIII B II 
Finally, analytic solutions of the linearized problem are found in Sec. IIII CI and studied in detail in Sec. IIII Dl Section 
IIVI contains concluding remarks and prospects for future work. The paper has two appendices, one of which contains a 
consistency check for the calculations performed, and the other describes in detail transition to the case of vanishingly 
small anchor dimensions within the large-slope expansion. 



II. PRELIMINARIES 



A. Spontaneous flame dynamics on-shell 



Consider a 2D-flame propagating in a channel of constant width 6, filled with an initially quiescent uniform ideal 
gas. Let the Cartesian coordinates {x,y) be chosen so that the channel walls are at x — 0,5, and y — — oo is in the 
fresh gas. These coordinates will be measured in units of the channel width^, while fluid velocity, v — {w, u), in units 
of the velocity of a plane flame front relative to the fresh gas, Uf. Finally, the fluid density will be normalized by the 
fresh gas density, 9 > 1 denoting its ratio to that of burnt gas. We assume that the flame pattern is continued to the 
whole X-axis in the usual way using the ideal boundary conditions at the channel walls: 

/' = 0, w = for x^O,b. (1) 

Then the on-shell value, (?/;_, tt_), of fresh-gas velocity (i.e., its value at the flame front considered as a gasdynamic 
discontinuity), and the flame front position, f{x,t), satisfy the following complex integro-differential equation [1, 0| 

f T 7 1 ' 

2 (u;^)' + (l + i^) S M - ^ / dS:{idy - d^) / drAf (.•?, i - t) I = , (2) 

I oo T- ) 

where uj = u + iw is the complex gas velocity, [to] its jump across the flame front, dx = d/dx, dy = d/dy, and the 
prime denotes differentiation with respect to a; (in the last term on the left, the argument y is understood to be set 
equal to f{x,t) after partial spatial difl'erentiation, but before the a;-differentiation denoted by the prime; we recall 
that the improper i-integral in this term is understood as an analytic continuation of the corresponding regularized 
expression, see Ref. for details). The memory kernel M has the form M{x,t) = N{x,t)v^{x,t)cr^{x,t), where 

— \J 1 + (/')^ J ^+ = ^i+n-i is the normal burnt gas velocity relative to the flame front, rii denoting the unit vector 
normal to the front (n points towards the burnt gas), and a-^- is the on-shell value of vorticity produced by the curved 
front. The memory kernel is integrated over any path in the complex time-plane, connecting the points 



where 



{w+,u+), u+{x,t) ^ u+{x,t) 



df{x,t) 



dt 

is the on-shell burnt gas velocity relative to the front, and r is the radius- vector drawn from the point [x, f{x,t)) at 
the front to the observation point [x, y). Finally, the action of the operator 3^ on an arbitrary fmiction a{x) is defined 
by 



(Oia) (x) 



+ 00 

l + z/'(.T,i) [ a{i) 

x^x + i[f{i,t)- f{x,t)]' 



^ However, we keep track of b throughout Sec. [IT] 
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where the slash denotes the principal value of the integral. For a 2&-periodic function a(x) [i.e., a{x + 2b) = a{x) 
summing explicitly the integrand with the help of the formula 



+00 

2bk + z " 26^°H26J 



/e— — 00 



the right hand side of Q can be rewritten as an integral over the channel width 



(j^a) {x 



l + if{x,t) 
2b 



dx 



(5) cot { ^(i - a; + t) - f{x, t)])} 



(4) 



We recall also that the value of vorticity at the front and the normal velocity of the burnt gas, entering the function 
M{x,t — r), as well as the velocity jumps at the front, are all known functionals of on-shcU fresh gas velocity [1, [^. 
For zero-thickness flame fronts one has 
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where 



Together with the evolution equation 



D _ d 



+ + 
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Dt ^ N Dt 



(5) 
(6) 



N dx 



• n) = 1 , (7) 
the complex Eq. Q constitutes a closed system of three equations for the three functions ix,t), u^{x, t) and f{x,t). 



B. On-shell description of anchored flames 

As derived, Eq. ^ describes only spontaneous flame evolutions. However, the anchoring system is not difficult to 
incorporate into the framework of the on-shell description. This can be done as follows. Consider the simplest and 
most commonly used in practice type of the anchoring system - a metal rod placed somewhere within the channel. 
From the mathematical point of view, the presence of the rod can be described as a singularity of the complex velocity, 
uj = u + iw, considered as an analytical function of the complex variable z = x + iy. Namely, suppose that the original 
field, u!o{z), is superimposed with the complex velocity, uj'^{z), describing a dipole located at the point {xo,yo): 

uJo{z) + '^—^ = uj{z) , (8) 

[Z — Zq) 

where zq — XQ+iyo, and d = di -\-id2 is a complex constant determining strength of the dipole as well as its orientation. 
For sufficiently small \d\, perturbation of the main flow is noticeable only in a small vicinity of the dipole. Since loq{z) 
is analytical at z = Zq, one has 

u)o{z) = ujq{zq) + 0{\z - zo\), 
and hence, the complex velocity near the dipole can be written approximately as 

uj{z) ujoizo) + -, (9) 



{z - Zq) 



The form of the stream lines is given by 



Re i uJo(zo){z — Zq) '> = const . 

Z- Zq ' 
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, \ r \ di{x - xq) + d2iy - yo) , 

Uo{X - Xq) - Wo[y ~ yo) r^— ^ = COllst , 

(X - xoy + [y - yoV 

where uq, wq are the real and imaginary parts of uJo{zo). It is seen that if we choose di = uqR^, c?2 = —wqR^, with R 
arbitrary real, then the stream-line family contains a circle of radius R, centered at the point (xq, tjq). Thus, adding the 
term uj'^^z) = ujq{zq)R'^ /{z — zq)^ to the velocity field u!q{z) describes perturbation of the given flow by a cylindrical 
rod of radius R, centered at zq. To take into account non- uniformity of the main flow near the rod, and to describe 
more general rod profiles, it will be necessary to superpose several dipoles located within the rod area, and to include 
higher-order multipoles into consideration. 

To obtain generalization of Eq. ([2]) to the case of anchored flames, we recall that this equation is a consequence of 
the following relations: 



(l-iOi) (w_)' = 0, (10) 
(i + zM)K)' = 0, (11) 



-t-oo 



ojI^- I dx{idy-dx) I dTM{S:,t-T) (12) 



-OO 



(13) 

where [lu] denotes the jump of the complex velocity across the flame front, [uj] = uj{x, f{x,t) -f 0) — uj{x, f{x,t) — 
0). Equations (|10p . pT|) express analyticity and boundedness of the complex velocity upstream, and its potential 
component downstream 0, Q , Eq. is the on-shell expression of the rotational component 7] , while Eq. is an 
obvious identity. As we have just seen, the presence of the rod violates analyticity of the complex velocity, so that 
either of Eqs. pO]) . (fTTj) is no longer valid, depending on whether the rod is placed up- or downstream. In the former 
case, Eq. ([T0| is satisfied by ujo{z) = uj{z) — u>'''{z), because it is analytical upstream and bounded. On the other 
hand, since uj'^{z) does not have singularities downstream and is bounded there, it satisfies Eq. (fTTj) . Thus, 



(l -^^C) {iu_-ujt)' = 0, 

(i + i^c) K)' = o. 



Since ojI = ojf^, we see that Eq. (flUl) is replaced in this case by 



l-i-K^ {lo^)' ^2{Loi) . (14) 
Accordingly, acting on Eq. by the operator (1 + i'K), we obtain the following equation 

r . y V 

2(w_)'+ (l + i3fW M- ^ / dx{idy-d^) j dTM{S:,t-T) \ = 2 {cut)' , (15) 

I. oo T- ) 

which is the sought extension of Eq. ([2]) to the case of anchored flames. In the case of the rod located downstream, 
similar considerations show that Eq. pT]) must be replaced by the following 

(l ~ i-Sdj [uj^)' ^ G . (16) 
It is not difficult to verify that the resulting equation for aj_ in this case has exactly the same form (1151) . 



C. Influence of anchoring system on V-flame structure 

As was mentioned in introduction, the necessity of explicit inclusion of the anchoring system into consideration raises 
the question as to what extent this system affects global properties of V-flames. Let us now show that as long as linear 
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dimensions of the rod are small compared to the channel width, so that the flame front can be considered piecewise 
linear, influence of the rod on the flame structure is local, in the sense that it is confined to a small region near the 
rod. We recall, first of all, that the relative value of the velocity disturbance caused by a dipole modeling the rod is 
proportional to B? /{x^ +2/^) (for simplicity, the dipole is assumed to be at the origin). Hence, under the assumption 
i? <C fo, this disturbance is indeed negligible for the most part of the channel, except a small region {x,y ^ R) near 
the rod. This simple reasoning is not yet sufhcicnt to prove our statement, because it only demonstrates the locality 
of, so to speak, direct rod influence on the flow structure. In such an essentially nonlocal problem as deflagration, we 
also have to look for possible indirect consequences of this influence, related to the fact that the presence of the rod 
ultimately determines the basic flame pattern. The on-shell description is particularly convenient for this purpose, as 
it explicitly reveals the nonlocal structure of the governing equations. 

For the rest of the paper, flames will be considered in the reference frame attached to the rod (the above-given 
formulation is invariant under transitions between different reference frames). Accordingly, the fresh-gas velocity at 
infinity will be denoted U: 

u{x, y = —oo, t) = U{> 0). 

We will assume in what follows that the anchoring system is stationary, i.e., its properties do not change with 
time. This means that these properties can be inferred from the steady-state V-flame structure. To this end, we note 
that the stationary version of Eq. reads (here we are in the rest frame of the flame- front, so the over-bar in the 
notation of velocity is omitted) 

2 (._)' + (l + {h' - = 2 (.1)' , (17) 

which follows directly from the fact that Eq. ([2]) reduces in this case to the stationary equation derived in [1, [E] ■ In 
regions where the flame-front slope is constant and the upstream flow is homogeneous, the first term on the left as well 
as the expression in the curly brackets vanish, because velocity jumps are constant there, uj- ~ const, and vorticity is 
not produced. This expression is only non-zero in a vicinity of the rod where all the quantities involved vary rapidly. 
It is this rapid variation that is a possible source of indirect influence of the rod on the global flame structure. Indeed, 
for R/b 0, both terms in the curly brackets have a ^-functional character. If the ^-singularity were not canceled in 
their sum, then upon the action of the !K-operator it would give rise to an expression which is non-zero everywhere in 
the channel. However, we have just seen that the right hand side of Eq. (|17p vanishes outside of small region around 
the rod. Therefore, in order that this equation be satisfied, the (5-contributions must cancel. To be more specific, let 
us assume that the rod is located downstream (which is normally the case in actual experiments) , as shown in Fig. [T] 
Then, using Eqs. ((TU|), (fTC)) . and the identity uj- + [uj] = uj+, Eq. p?)) can be conveniently rewritten as 

There are two types of (5-like contributions on the left hand side of this equation, corresponding to the real and 
imaginary parts of the expression in the curly brackets. Since the real part is even under x — > — x, its x-derivative is 
odd. Hence, the corresponding singularity is generally proportional to S'{x), and can be compensated by appropriately 
choosing the coefficient d in the dipole term on the right hand side. Indeed, the on-shell value of a dipole uj°'{z) = 
a/ [z — a)^, considered in the limit a — + 0, possesses all characteristic properties of the (5-function: lo+{0) — \/a~^ oo, 
uj+{x) — > 0, for X 0, and the integral dxa/{x + if{x) — a)^, taken over a region A 3> a around x = 0, has a finite 
value (because f{x) is an even function). 

Things are different, however, for the imaginary part which is odd in x. In this case, the singularity is proportional 
to 5{x); for zero-thickness flames, for instance, contribution of the first term in the curly brackets to the singularity 
is equal to — 2i(^? — l)s5{x) / \/l + s^, where s is the value of the front slope far from a; = 0, as is seen from Eq. ([5]). 
Singularities of this kind^ cannot be compensated by any local field w'^(z). Thus, we arrive at the conclusion that the 
assumption of piecewise linear front structure implies the absence of terms proportional to 5{x) on the left hand side 



^ In fact, it is the singularities ~ with undifferentiated (5-functions, which are only important. Indeed, on dimensional grounds, a 

differentiated 5 should be accompanied by an extra factor with the dimension of length; since this an "inner" contribution, the factor is 
~ R, and hence the 5' (a;)-terms can be neglected in comparison with &{x) in the limit R ^ 0. Another way to see this is to recall that 
the parameter a in the dipole u}'^{z) = a/(z — a)^ is 0{R), while the strength of the dipole modeling the rod, \d\ = O(R^), as we saw in 
Sec. Ill B1 Hence, uj'^{z) must be accompanied by a factor 0{R). 
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of Eq. (fT8|) . i.e., that the contribution of the first term in the curly brackets is canceled by that of the second term. 
This requirement can be written in the following integral form 

where A: i? <C A ^ 6 is the length scale where "inner" solutions (|a:| <^ b) are to be matched with the "outer" 
ones (|a;| ^ R). Indeed, by virtue of Eq. the contribution of the small region near the rod to the left hand side 
of Eq. is also small outside that region, which is just the required absence of the (5-terms. Equation thus 
represents a condition that selects inner solutions compatible with the prescribed global flame structure. 



= 0(1) for A/b~^0, 



(19) 



D. Linearized equation for flame perturbations 



In the present paper, we for are looking for possible genuine V-flame instabilities, which would be inherent to the 
V-configuration itself, and unrelated to the properties of a specific anchoring system. We thus assume, as was already 
mentioned, that this system is stationary, and the condition p^ is fulfilled. Then the equation for flame perturbation 
is obtained by linearizing Eq. psp around the stationary solution, with the right hand side kept fixed. This linearized 
equation thus coincides formally with that derived in Ref . 0, 01 , but for our present purposes another form of this 
equation will be more appropriate, which avoids explicit differentiation of the memory kernel. 

First of all, since the basic pattern is stationary, time-dependence of perturbations factorizes: 

Sf{x,t) = f{x)e''\ Sw-ix,t) ^ w{x)e''\ Su^{x,t) ^ u{x)e''\ (20) 

where v is a complex constant to be found as part of the solution. Not to mix the imaginary unit entering v with 
that appearing in Eq. IjlSp . we will denote the former by j: 

where U12 are real numbers. Accordingly, the amplitudes /, w, u are to be understood complex with respect to j 
(until Sec. IIII Di j will not appear in formulas explicitly; an example illustrating the use of this "double imaginary 
unit" formalism is given in Appendix A). Next, taking into account that the basic solution is piecewise constant, we 
obtain the following equation for the x-dependent parts of the perturbations 




i j diM„(i)|„(i)^exp(^-^e-*^jx(^-i) I =0, (21) 



where {^a} = {f,'w,u), Ma{x) is the differential operator obtained by linearizing the function M{x,t) around the 
stationary solution, and setting d/dt ^ v afterwards; (/) € [— tt, -|-7r] is the angle between the vectors r, defined 
positive if the rotation from to r is clockwise, and x{x^ is the sign function, 

-1, x > , 
Xix) ={ 0, a; = 0, 
-1, a; < 0. 

In the form written, Eq. (j2ip applies to flames with arbitrary jump conditions at the front and arbitrary local 
propagation law. However, to investigate the problem as stated in the beginning of this paragraph, we do not need 
to remain at such a general level. As mentioned earlier, we are concerned with instabilities specific to the presumed 
V-pattern, so the characteristic perturbation wavelength of interest is of the order of the channel width which is 
normally much larger than the cutoff wavelength. Hence, the curvature effects can be completely neglected in our 
investigation, and the consideration be limited to the case of zero-thickness flames. Then the linearized velocity jumps, 
appearing in Eq. (|2ip. take the form 

To linearize the memory kernel, another form of the expression ([S]) will be more suitable, which avoids appearance of 
the second spatial derivatives of the flame-front position. The point is that linearizing Eq. ([S]) directly is readily seen 
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to lead to expressions of the type x{x)5{x) which are not well-defined in the sense of distributions. To resolve this 
ambiguity, we first rewrite Eq. ([7]) as 

u--%- f'w- = N , 
dt 

differentiate it with respect to x, and use the resulting equation to eliminate /" from the right hand side of Eq. ([6]). 
The memory kernel thus becomes 



M = -(6-1) 



J 1" W-W_ + U-U_ 



dt ' dt - dt 



(23) 



The right hand side of this expression now involves only first derivatives of continuous functions. It should be 
emphasized that this trick works only in the outer region where effects due to the finite flame-front thickness are 
negligible. If they are not, /" appears already in the undifferentiated evolution equation. Linearization of Eq. (j23[) 
yields 



Maic{x) = -{9 - 1) iyw{x) + sx{x)vuix) + vT+7^u'{x) . (24) 
Finally, the linearized evolution equation reads 

u{x) - sx{x)w{x) = v~f{x) + f^M£M . (25) 

Vl + s 

1. Boundary conditions 

We consider symmetrical basic V-patterns, so that the flame-anchoring rod is located in the middle of the channel. 
For simplicity, flame disturbances also will be assumed symmetrical under reflection with respect to the y-axis. In 
these circumstances, it is convenient, without changing notation, to consider a double-width channel occupying the 
strip —h X ^ with the rod being at the origin a; = y = 0, and the y-axis playing the role of the symmetry 
plane of the flame. Accordingly, boundary conditions of the exact problem, i.e., the problem we started with, read 
w = f = for X = ±&. They are used, in particular, for the periodic continuation of the flame pattern, mentioned 
in Sec. Ill Al After the initial problem is divided into the inner and outer ones, these conditions naturally remain 
pertaining to the former. A question thus arises as to the boundary conditions relevant to the outer problem. 

Evidently, the requirement of a vanishing front slope at the walls is now irrelevant. Indeed, it is not satisfied already 
by the basic steady V-configuration defined to have the constant slope, |/'| = s, everywhere. The front flattening takes 
place in a thin boundary layer near the walls, characterized by large values of /". On the other hand, the a;-derivative 
of the ui-componcnt does not have to be large in this region, as can be seen from the fact that the condition w = 0, 
being the universal boundary condition for the ideal fluid, is compatible with any prescribed front configuration. 
Indeed, the linearized Euler equations for the fresh-gas read 

^ ^dSw d6p 

+ = (26) 

J. ^dSu dSp 

where it is taken into account that for the steady-state V-fiame, u = U, w = Q. Using the continuity equation in 
Eq. (I27|) . multiplying it by /', subtracting from Eq. (j26p . and going over on shell one obtains an equation for 5w^: 



or 



dSw- 1 I / c- N f d5p\ 

-ir = uY^'^-^^^[-d^)_ 



(28) 



where I is the front length counted off from the channel wall, and n is, as usual, the unit vector normal to the front. 
As was already mentioned, the front flattens in a thin boundary layer near the walls. From the standpoint of the 
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outer problem we are concerned with, a finite change of an on-shell variable across the layer is seen as a finite jump 
of that quantity at the channel wall. Accordingly, its derivative with respect to / contains a term proportional to 5{l). 
Now suppose that this is the case for 6w-{l). Then it follows from Eq. (|28p that 

where c is a constant, and "• • • " denote regular terms. Since the velocity jump is finite, this means that 

This relation can be integrated by noting that differentiation of pressure in the direction normal to the front cannot 
produce a (5-singularity along the front, and hence, 

Sp ^ 6(5(0 + ■■■ , 

where 6 is such that (9C/9n)_ — c. But pressure is only allowed to have a finite jump at the wall, therefore, 6 must 
vanish upstream. Meanwhile, in the absence of obstacles at the wall, the outer solution is regular on-shell, and hence 
the function C is differentiable not only in the near upstream region (as required by its definition), but also at the 
flame front. Under such circumstances, the requirement 6 = upstream entails vanishing of its derivative at the 
front, i.e., = 0. Thus, Sw- is in fact continuous at the wall, so its vanishing remains a boundary condition of the 
outer problem: Sw^ (±6) — 0, or 

w{±b) = 0. (29) 

The reasoning just given does not apply at cc = 0, because of the presence of the rod. Nevertheless, w{0) must also 
vanish, as a consequence of our assumption that the anchoring system is stationary. To see this, let us consider the 
procedure of matching the inner and outer solutions in more detail. Take the x-component of the fresh-gas velocity. 
For gas elements moving near the y-axis, this component is zero everywhere except a small vicinity of the rod. More 
precisely, w induced by the rod is 0{U) for \/ + -tp = p ^ and rapidly decreases with distance. At distances 
p ^ Rq, where R <ti Rq ^ b, the inner solution describing the flow near the rod is matched with the outer solution 
we are interested in. In the steady case, matching at the flame front assigns W- a definite value, say wq, which is 
generally nonzero. This value plays the role of a boundary condition for the steady flow, deflning thereby the basic 
pattern. Now, since the properties of the rod are assumed stationary, in particular, unaffected by perturbations of 
the outer solution, so is the flow near the rod. Therefore, matching of the inner solution with the outer one will give 
w the same value Wq. In other words, 5w-\pr^jig — Swq = , which in the limit i?, i?o ^ yields 

?i)(0) = 0. (30) 

By the same reasoning, 

w(0) = 0. (31) 

Finally, the remaining condition replacing /' = is 

/(0) = 0. (32) 

It follows directly from the fact that we consider the rod dimension as vanishingly small compared to the channel 
width. Indeed, the linear dimension of the fiame tip as well as its separation from the rod are both ~ R. Hence, 
f{x ~ i?) ~ i?, which in the limit R ^ gives Eq. (|32p . This condition means that the flame is not torn off from the 
rod by the perturbation. 

III. ON-SHELL DYNAMICS OF V-FLAME PERTURBATIONS 
A. Evaluation of the rotational contribution 

In order to study evolution of the V-fiame disturbances using Eq. (I^T]) , we have to evaluate the improper x- integral 
appearing in the curly brackets. We recall that this integral is understood as an analytic continuation of the regularized 
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integral 

+ CXD 

J die-^l*-^lif„(i)C„(i)^exp X{x - i) , (33) 

—00 

to the limit /i ^ 0^. To simplify the calculation, we note that 

rw+e^*"^ = — i(z — z)liJ^ , z = x + is\x\ . 

Indeed, one has \z — z\ = r, \uj^\ = v^, while according to the definition of the angle (f> it is equal to the phase difference 
of the complex functions w+ + iu+ = iu^ and (z — z). We need to consider two different situations corresponding to 
the integration variable running over the negative- or positive-slope part of the flame front (see Fig. [2]). Assuming 
that the observation point x E [0, 1], one has, in the case x £ [— 2n, —2n + 1], n £ Z, 

z — z = (x — 1] + 2n) + is{x — rf) , 
where [0, 1] 3 i] = x + 2n . Similarly, in the case i G [—1 — 2n, — 2n], 

z — z = {x — rj + 2n) + is{x + rf) , 
where [—1,0] 3 r] — x + 2n . In effect, the exponent in the integrand of ([55)1 takes the form 

[~i{x -T] + 2n) + s{x-i])]\ , 77 e [0, +1] , 
[-i{x - Tj + 2n) + s{x + ri)]) , tj e [-1,0], 



exp e 




1^0 1 
where 

1 



Furthermore, the regularizing factor e ^1 may be replaced by e ^^'"1 , because (x — 77) is finite. Next, the i-integral 
taken over (— cxd, +00) can be represented as an integral over 77 G [—1, +1] of the integrand summed over all n. Since 
the functions Ma{x), uj+{x) are periodic by construction, we need to sum the following series 



+00 



/(/i) = exp {27T,ix — 2|n|/i} x(a; — ?7 -l- 2n) , 

where 



n— — 00 



v/ujo, Tje [0,+l], 
iy/Lo*o, [-1,0]. 



Taking into account that [a; — 77I ^ 2, one has 



+00 -t-00 
I{fi) = x{x — ^?) + 6xp {2n{i>c — fi)} — exp {2?i(— — /i)} 

n=0 Ji=0 

= Xix - ?7) + [1 - cxp{2(z>f - ^i)}]-' - [1 - exp{2i-i^ - fx)}]-' . (34) 

Since the initial improper i-integral is reduced to an integral over a finite domain, its analytic continuation to = 0+ 
amounts to that of the function /(/i), which is 

/(0+) = — ^7) + [1 ~ exp{2i>f}] "'^ — [1 — exp{— 2i>f}] ^ — x{x — rj) + i cot k . 

All these formulas were derived for x £ [0, +1]. From these, the corresponding expressions for x G [—1, 0] can readily 
be obtained by noting that the integral (|33p is invariant under the combined operations of inversion x —x, and 
complex conjugation. This rule can be deduced directly from the explicit formulas (|22p . (|24p . taking the various parity 
properties of the flow variables into account, yet it is in fact a general property of the formalism, unrelated to the 
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specific approximations made. In what follows, we will denote this combined operation as {x — > —x)*. It should be 
kept in mind that the complex conjugation here is understood with respect to the imaginary unit i, but not to j: 

■ ■ * 

Putting all these results into Eq. (j2ip , we thus arrive at the following linearized equation governing evolution of the 
flame disturbances 



2(D' + 



- 1 



^i>i:(x-\-is\x\) 



UJQ 



+ 1 



dr] vw{r}) + svu{rii) + \fl~\~^u {-q) 



xe 



(35) 



where the symbol (x — > —x)* refers to the whole expression written out explicitly in the curly brackets. As a useful 
check of the calculations performed, it is verified in Appendix A that in the particular case s — Q this equation 
reproduces the well-known Darrieus-Landau dispersion relation [ll| for the perturbation growth rate. 



B. The high- velocity limit 

In its general form, Eq. (|35p can presumably be solved only numerically. It turns out, however, that it is amenable 
to a full theoretical analysis in the case when the velocity of the incoming fresh-gas flow is high: 

C/> 1. 

Being opposite to that of classical analysis [l^, [ll|, [l^, [l3| , this limit is of considerable interest both from practical 
and theoretical points of view, as it represents the situation where propagation of the flame disturbances is strongly 
affected by the basic flow. We will see that the nonlocal interaction of flame perturbations with the background takes 
a new form which is principally different from that encountered in the conventional weak-nonlinearity analysis. Also, 
dependence of solutions on the gas expansion coefficient becomes quite intricate, having nothing in common with that 
found in the small-gas-expansion approximation. 



1. Large-s expansion of the M -operator 



We start discussion of the high-velocity limit by deriving an approximate expression for the 3^-operator appearing 
in Eq. ([55)) . There, it is defined at the unperturbed front, f{x) — s\x\, 

+1 

(O^a) (2;) = i^t^l^ ^dia(i)cot{|[i-a; + is(|i|-|x|)]} . (36) 
By virtue of the relation 



large values of U imply that the front slope is also large, so the argument of cotangent in Eq. ([36| has a large imaginary 
part for almost all values of the integration variable, in which case one has 

cot {| [i - x + isi\x\ - |x|)]} « -ix{\i\ - \x\) . (37) 

This approximation is valid for all x except two small regions near x — ±\x\ . More precisely, taking into account that, 
for real ai,2, 

g(o2— mi) _|_ g — (02 — iai) , , 

cot(ai + ^2) = ~i —, = -ix(a2) + O ( e^^l"^' ) , 

^ ' g(a2-Mi) _ g-(a2-iai) \^ y ' 
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we see that Eq. (j37l) holds true, with an exponential accuracy, everywhere except 

5 : 1^1 e (|a;| - S, \x\ + S), 

where S — 0(1/ s). 

To develop an asymptotic expansion of M in powers of 1/s for s ^ 1, let us choose a real £ > satisfying 

e < 1 , s£ > 1 . 
Then the integral in Eq. (|36p can be rewritten as 



(38) 



a(x) cot |— [x — X + — , 



X — € X — € +1' 

+ 

— 1 —x-\-6 x-\-e 



dx a(x)x(|x| — x) 



-x-\-e x-\-e 



—x—e x—e 



(39) 



where we assumed that a; > 0, for definiteness. Notice that in the last term on the right hand side of Eq. ((39)l . only 
one of the two integrals is defined in the principal value sense. As such, it is proportional to the derivative of a{x). 
It is not difhcult to verify that contributions of this kind give rise to terms of the order 1/s^. Below, we will need IK 
expanded up to 0(l)-terms, so the principal-sense integral can be neglected. The other integral can be evaluated as 
follows, within this accuracy, 

-x+e +£ 

j dx a{x) coi^—[x — X ~\- is[\x\ — x)]^ — —ia{~x) j dx coth |— i + Trixj- 



+7rs£/2+7rix 



'ia{—x)- 



dy coth y . 



(40) 



-7Tse/2+TTix 



By virtue of the conditions ([38]), the y-integral can be calculated, with exponential accuracy, using the contour 
deformation shown in Fig. [3] 



+7rs£/2+7r't. 



dy coth y 



-■Kse/2-\-7rix 



— 7rse/2 -\-7Tse/2-\-7Tix 



+ TTSS/2 



J + J dy + -j- dy coihy ~ in — 'Ki{2x — \) 



-■Kse/2 



-■Kse/2+7rix +-Kse/2 

On the other hand, replacing cotangent by the sign function gives zero within the same accuracy 

—x+e +e 

= 0. 



J dx a{x)x{\S:\ — x) — a{—x) J dx x(i) 



Using these results in Eq. ([39| . and then substituting it in Eq. ([36]) gives finally 

+1 _ _ 

(Wa) (x) = (sxix) - *) fdi «(^)+^«(-^) ^(^ _ 1^1) + ,a{-x){2\x\ " 1) + O ( J 



(41) 



where the symmetry of the operator i3i under {x —x)* was taken into account to dismiss the condition x > 0. 
As a special case of this formula, let us consider the action of 3i on a derivative. If a{x) satisfies a(0+) = a(O^), 
a{+l) — a(— 1), then integrating by parts in Eq. (|4T|) readily gives 



Ka' ) {x) ^ {sx{x) - ^) {-a(|.T|) + a(-|.T|)} + ta'{-x){2\x\ - 1) + O ( - 



(42) 
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where the prime now denotes the derivative of the function with respect to its argument, a'{y) = da{y)/dy. It turns 
out that this formula holds true even if the function a{x) does not satisfy the above conditions of periodicity and 
continuity at the origin. This is proved in Appendix B. 

To conclude this section, some comments concerning the structure of the expression (j4ip are in order. First of all, 
it is seen that the result of the action of depends essentially on parity properties of the function a{x), namely, 
'Ka — 0{s), if a{x) is even, and !Ka — 0(1), if it is odd. Next, the appearance of a term proportional to a{—x) 
encodes a peculiar interaction between the points x and — x, which is natural taking into account that the front wings 
get close to each other in the limit s ^ oo. Finally, it should be noted that although the identity IK o IK = —1 is valid 
whatever the shape of the flame-front, in particular, in the large-s limit, it cannot be verified using the expression 
on the right of Eq. (I4ip . already because of the composition of its leading term with the undetermined remainder 
0(s)oO(l/s) = 0(1). 



2. Equation for the x-component of velocity. Relative order of the flow perturbations 

The results of the previous section allow us to a obtain a simple relation between components of the perturbed 
on-shell velocity. We use the following equation 



1 - i^j Cb' ^ , (43) 

which is obtained acting by (1 — iJ-C) on Eq. ([35ll . and using the identity IKo Jf = —1; it can be derived also directly 
by linearizing Eq. PH)) . Applying the formula ([1^ yields 

^^—^u'ix). (44) 
s 



We see that the boundary condition (|29|l is met explicitly, while setting x = and using ([30]) leads to a new condition 

m'(0) = 0. (45) 

It will be shown in the next section that this condition is also satisfied automatically by the solutions of Eq. (|35p . 

Next, we use Eq. (j44p to determine the relative order of the flow perturbations within the large-s expansion. It is 
convenient to assume that u — 0(1). It follows then from Eq. that w — 0(l/s), while using these in the linearized 
evolution equation ([25]) tells us that / = 0(1). Applying these estimates to Eq. ((35)) shows immediately that the term 
{i + sx)f'/{^ + s^)^^^ in the curly brackets can be omitted. Since the 7y-integral is explicitly continuous at a: = 0, so 
is the expression in the curly brackets, as was to be shown. 

In connection with this result, it is worth mentioning that the term (i + sx)f' /{I + s^)"^/^ represents the linearized 
velocity jumps which define the Frankel potential- flow equation That this contribution is negligible means the 
evolution of disturbances in the case under consideration is essentially rotational, and cannot be described within the 
potential-flow model. 

C. Analytical solution of the linearized equation in the high- velocity limit 

We are now in position to proceed to analytical solving of Eq. ([55]) in the case of high stream- velocity. Although the 
following calculation is a straightforward application of the formulas derived in the preceding section, it is somewhat 
lengthy. We give it in considerable detail because some of its points are definitely worth to be mentioned. 

1. Derivation of the integro-differential equation 



To begin with, it is convenient to rewrite Eq. ([35]) as 



2Cj' +^-^[l + i-K]E' = (46) 



E{x) = / drj v'w{rf) -\- svu{ri) + yl + ^u' [rf) 







X e-*-(i+*^)'' [i colx + xix- ry)] + (x ^ -x)* . (47) 
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The order of the leading contribution to the left hand side of Eq. ([46)) can be read off from its first term, uj' . According 
to the estimates of the previous section, it is 0(1), and is contained in the real part of the equation. To extract 
the relevant contribution from the integral term, we recall that the action of 3f on odd and even functions gives 
rise to terms of the order 0(s) and 0(1), respectively. Furthermore, taking into account that wq = 0(s), and hence 
K = 0(l/s), one sees that E{x) = 0{s). Therefore, according to the naive power counting the integral term is formally 
O(s^). However, there is actually no discrepancy in the orders of the two terms, because the 0(s)-contribution to 
E{x) turns out to be imaginary ewen, and thus cancels with its counterpart from (x — > —x)*. Yet, the formal estimate 
means that expanding imaginary part of E{x), one must generally keep terms up to the second relative order in 1/s. 
With this in mind, we write 



1 



1/2' 



(48) 



and then 



^i>c{x — l)—>cs{\x\ — l) 



(x-1) 



15 + 



iiy{9- 1) 



1/2), 



(N-1) 



\x\ - 1) 



(49) 



On the other hand, since in the factor [vw + svu + vl + s^{t'] all terms are real, it can be replaced by s{vu + u'), 
with no risk of mixing orders. Similarly, one can replace loq^ coi{v /ljo) by 1/z/, because the imaginary correction is 
0{l/s'^). Also, before expanding, it is convenient to integrate by parts the term proportional to u'{rj). Taking into 
account the boundary condition pip , we thus find 



E{x) 



su{l) 











)+- 


u{ 




/ t^o 




i 

- + 


X{x - 


v) 


V 







where 9{x) is the step function, 



9{x) 



+ 1, a; > , 
0, x^O. 



Expanding further within the required accuracy with the help of Eqs. (j48p . (j49p . and omitting contributions which 
are real odd or imaginary even gives 



E{x) = e''(i-l^l)M(l) 
1 



a{l~\x\) + 



ixia — V) 



2u{x)0{x) 1 



la 
s 



-(a + 1)6-^^1^1 / d??7i(??)e'^" ( 1 + ^^"^ "^^^ ) + (^ 



where a = 61 — 1. It is seen that the odd contributions arc of the order 0(l/s) indeed, so upon the action of !K they 
give rise to 0(l)-terms. 

Extracting the real part of Eq. (|46)) with the help of the formula (142]) gives 

2u'{x) + (6* - l)\x\KeE'{x) + s{e - l)x{x)lmE{\x\) = . (50) 

Since E{x) is given by an integral of a piecewise continuous function [Cf. Eq. (j47|) ]. it is continuous. Therefore, its 
imaginary part being an odd function turns into zero at the origin. Then Eq. IjSOp tells us that its solutions satisfy 
the boundary condition flS)) . 

Substituting the above expression for E{x) in Eq. ()50p . and introducing a new unknown function g{x) according to 



u{x) = g(x)e 



-v\x\ 



(51) 
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we finally obtain the following integro-differential equation 



\x\ 

(1 + a\x\) g'{x) — {o? + OLv\x\ + v)g{x)x{'x) + a(a + \)vx{x) J drig{ri) 



+aiygix{a\x\ — a — 1) = , 



where gi = g{l), and we used the identity 



1 

dV9{v)[x{\x\ - r/) + 1] = 2 J df]g{r]) . 





(52) 



2. Solution of the integro-differential equation 



Up to an additive constant, Eq. (|52p is equivalent to the following ordinary differential equation obtained by 
differentiation with respect to x [in view of the symmetry of this equation under a; — s- —x, it is sufficient to consider 
it on the interval x G (0, 1)] 

(1 + ax) g" {x) — ic? + avx + v — a)g' (x) + a^i>g{x) + ai'gi{2ax — a — 1) = . (53) 

The general solution of this equation can be found in the form 

g{x) ~ Ci + C2VX + h{x) , (54) 

where Ci_2 are constants, and h{x) satisfies 

(1 + ax) h"{x) - (a^ + avx + v- a)h'{x) + a^vh{x) = . (55) 

The latter equation can be reduced to the degenerate hypergeometric equation, and its general solution conveniently 
written as 



h{x) = 03 2; + 



x+l/c 



dyy 



(56) 



fS/a 



where C3 and /3 are new constants. A direct substitution shows that (j54p is a solution of Eq. ((53|) . provided that the 
constants P, Ck, k — 1,...,3 satisfy 



(a^ — a + y)i'C2 — a^uci + a{a + Vjvgi = , 
[a — l)vc2 + 'iagi — . 



(57) 
(58) 



In addition to that, for ([M]) to be a solution of the integro-differential equation ([5^ . the constants must be chosen so 
as to guarantee vanishing of the additive constant in this equation, which was lost upon the transition to Eq. (j53|) . 
To extract this constant, we first of all note that 



(1 + ax) h'ix) - ia^ + a.x + .)/.(.) + a(a + 1). / d.H,) = 



V j dyy 



which can be checked by direct computation. Then collecting the additive constants in Eq. (j52|) gives another equation 
for P, Ck- 



VC2 - (,a + vjci + C3 



ae 



u I a 



V f dyy e 



(59) 



Analytical stability analysis of V-Bames, v 1 - February 26, 2009 



15 



Finally, the boundary condition ((3T|) takes the form 



C1 + C3 J dyy-'^-^e'^y/" = . 



(60) 



Four equations ([57)1 - ([50)1 constitute a closed system for the four constants /3, Cfc. In particular, the condition of 
consistency of this system determines the spectrum of the perturbation growth rate v. The boundary value of 5(2;), 
entering these equations, is expressed through the unknowns as 



gi= ci+ C2V + C3 / dyy 

/3/(a+i; 



-CK — lgi/(a+l)y/a 



(61) 



3. Reduction to an algebraic system of linear equations 

Since Eqs. (I57p - (j60p were derived from relations linear with respect to g{x), by an appropriate redefinition of 
the unknowns they can be naturally rewritten as a system of linear homogeneous equations. For this purpose, let us 
introduce the following notation 



It is not difficult to check that 



1 

/ 

/3 



dyy 



V 

n = - , 



C4 = C3$[n,/3], 
1 



$ = 



dyy 



-Q-lg(a+l)ray 



(62) 
(63) 
(64) 



l/(a+l) 



/3/(q + 1) 

Using this in Eqs. (j60p . (|6ip allows us to put them into the form that no longer involves /3 explicitly: 

ci + C4 = , gi = ci + C2i^ + C3<1) + C4(a + 1)" . 

On the other hand, since gi is linear with respect to Ck, k — 1, ...,4, so are Eqs. ([57)1 ~ ([50]) . Therefore, taking C4 as 
an independent unknown instead of /3 renders the system linear algebraic. Eliminating C4, we thus obtain 



Cia + C2 |(a — 1) 
C2(a - 1) 



a + 1 



2 

a + 1 



n - 1 



-ri^ = 0. 



n- 1 



ci [1 (a + 1)"] +C2(3a- 1)- +C3$ = 0. 



(65) 



D. Structure of the solution 



1. The perturbation growth rate spectrum 

The solvability condition for the system reads 

<I>e""a(a - 1) [{a + l)n - 2] - n {a{3a - 1) + [{a + 1)" - 1] {a^ - 3)} 
+2(a-l)[(a + l)"-l] =0. 



(66) 
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This equation determines the spectrum of flame disturbances, i.e., the admissible values of the perturbation growth 
rate, v. Before looking for its numerical solutions, it is useful to establish general features of the spectrum. For this 
purpose, it is convenient to switch from a back to 6 — a + so that the definition (|64p takes a more compact form 



/ 



$ = \ dyy 



eBny 



Integrating by parts, we can rewrite this formula for |n| ^ 1 as 

$ = -1 {e«" - fl^'e"} [1 + 0{l/\n\)] , \n\ > 1. 
On 

It is evident from this expression that <I> ~ e^" for Re n +oo, and hence has no solutions for such n's. On the 
other hand, <!> ~ e" for Re n ^ — oo, which is compensated by the factor e~" in Eq. ((66)) . However, the coefflcient of 
the combination $e^" as well as the rest of the equation are polynomials in n, so there are no solutions in this domain 
either. Thus, eigenvalues tend to be vertically aligned in the complex plane. Substituting the above asymptotic into 
Eq. (Uni) yields 

fl (6* -l)(36i- 4)+ [6*^-1 -ll (612-261-2) 
e'^ = 6" + S{9),. , S{e)^'- ^ (0 - 1)2(^0 _ 2) 1^1 (^7) 

Despite appearance, the function S{d) has no pole at 6 — 2 (see Fig.Hl). 

As we just mentioned, the simplified relation (j67p determines the spectrum in the case |Im iy\ 3> 1. From the 
practical point of view, however, we are interested in v^s whose imaginary part is not too large, so that only a finite 
number of eigenvalues need to be taken into account. Indeed, recalling the relation u(x) — g{x)e~^'^ , the characteristic 
wavelength of flame perturbation with the given v is 

27r 



Im V 

In terms of displacements along the front, AZ — sAx, this corresponds to a wavelength 

27rs 



A 



Im V 



On the other hand, perturbations with wavelengths less than the cutoff wavelength, Ac, are damped by the curvature 
effects. The condition A > Ac gives, in ordinary units, 

^ 27rte , , 

Im 1/ < — — . (68) 
Ac 

For gas expansion coefficients of practical importance [6 = 5-;- 8), the quantity 9^ is very large; S{9) is also large, but 
smaller than 9^ by about two orders. It follows from Eq. ((67)) that if imaginary parts of the eigenvalues are not too 
large, they are close to multiples of 27r, while their real parts are approximately equal to 9\n9, 

i^m = 9lii9 + j2Trm, meZ, 6* » 1 . (69) 

This formula is useful for searching and identifying numerical solutions of the exact relation ()66p even for smaller 
values of 9. Its validity as a classification scheme breaks when S{9) « 9^ . In fact, purely real solutions exist for 
9 < 9q Ri 1.8. The corresponding modes describe aperiodic development of disturbances. 

Examples of I'-spectra obtained by solving Eq. ()66p numerically are presented in Table I. Figure [5] illustrates 
graphical determination of the lower parts of n-spectra. They show that all solutions have positive real parts. 

To conclude, for sufficiently large values of the incoming fresh-gas velocity, the piecewise linear V-structure is 
unstable for all values of the gas expansion coefficient. 
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A 
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iD.ol + Jll.oy 


J IZ.OD 


n 
6 


/111 1 ^'on on 
4.ii + J Z\j.Z\j 


y.i / + J 19. zo 


Id. <2 + Jlo.21 


Jlo.oO 


4 


4.37 + j26.52 


9.47 + J25.76 


17.05 + j24.85 


J25.13 


5 


4.58 + j32.83 


9.70 + j32.20 


17.32 +j31. 39 


i31.42 


6 


4.76 + j39.14 


9.89 + j38.60 


17.54 + j37.86 


i37.70 


7 


4.90 + j45.44 


10.05 + j44.96 


17.73 + j44.29 


j43.98 


8 


5.03 + j51.73 


10.18+ j51.31 


17.89+ j50.70 


i50.27 


9 


5.15+ j58.02 


10.30+ j57.65 


18.03+ j57.08 


i56.55 


10 


5.23 + j64.31 


10.41 +j63.97 


18.15+ j63.45 


j62.83 


ein6 




9.4 


18.2 





TABLE I: Lower parts of the perturbation growth rate spectra obtained by solving Eq. (I66p numerically. The eigenvalues are 
measured in units Uf/b. The last row and the last column list their real and imaginary parts as given by the formula (|69|l . 



8. Space-time profiles of the flow perturbations 

To write down solutions for the flame perturbations, we need to represent Eq. (j54p in a form suitable for separating 
its real part. Using the definitions ([62)) . ((63l) . one has 



l+ax 

a;+l/ a 

dyy' 

l/a 



1 1+ax 



dyy 



e " 



C4 (1 + aa;)" + C3 H ^ y 



(70) 



Also, the complex phase of one of the coefficients Ck in the linear problem can be chosen arbitrary. We use this to 
make C3 real. Then, writing i' = +ji'2, Ck = \ck\e^'^'' , one combines the formulas (|20p . (|5ip . (IM)) . (|56|) . and extracts 
real parts of the resulting expressions. Thus, we find 

5u^{x,t) = \ci\ [1- (l + ax)"]e''i(*-^)cos[z/2(<-.x) + (/:'i] 



\C2\xe 



{vi cos[i^2(t - x) + (^2] - ^2 sin[j/2(t - x) + (^2]} 

-1/q 



+C3 + ^) / d2/2/^"-ie'^i(«+*~^) cos[i/2(2/ + t - x)] , 



(71) 



l/a 

The corresponding expression for the w-component follows then from Eq. 

1 — X d 



5w_ [x, t) 



-5u^ {x, t) 



s dx 

Finally, in terms of the function .9(2;), the linearized evolution equation (|25p takes the form 



(72) 



[e^^fix))' = g{x)[l + - x)] + (x - l)g'{x) 
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Substituting the solution (|54|) . and integrating gives 



x + - 
a 



l + -)a;+--2-l 



a + 2 

+a;e~''^[ci(l + i^) - 02^] + x^e^'^^z^ 

where C5 is a constant. Its value is fixed by the condition p2l 

1/ a + 1 



V I C 



a + 2 
2 



1 1 
J/ a 



056 



C5 = Ci 

The perturbed front shape is given by 



a2 a + 2 



1 



e''/" /I 


1 




+ ^^^ + 2^^ 






a 





(5/(x,i) =Re{/(x)e''*} , 



(73) 



which we do not write out explicitly because of its complexity. 

All expressions above are written for a; > 0. They can be easily continued to a; < using parity properties of the 
flow variables. 



IV. DISCUSSION AND CONCLUSIONS 



The results of analytical investigation presented in this paper give an accurate and complete account of the stability 
properties of confined V-flames anchored in high-velocity streams. The general conclusion we arrived at is that in 
this case, the piecewise linear V-structure is unstable for all values of the gas expansion coefficient. The perturbation 
growth rate spectra have a similar structure for all ^, obeying simple classification with respect to the imaginary part 
of eigenvalues. The only exception is the existence of aperiodic unstable modes for flames with 6* < ~ 1-8. We have 
found also explicit analytic expressions for the eigenfunctions [Eqs. ([7T|) ~ ((73)) ]. 

One result that deserves special emphasis is that dynamics of flame disturbances in the high-velocity limit turned 
out to be governed by the memory effects associated with vorticity generated by the curved front, which completely 
dominate contributions due to gas- velocity jumps across the front that define flame behavior in potential models. This 
is in striking contrast with what has been found for freely propagating flames, where development of the Darrieus- 
Landau instability is determined mainly by the structure of these jumps (in the Sivashinsky-Clavin [Tsj and Frankel 
[13] models, for instance, memory effects are completely neglected). 

Furthermore, dependence of the solution on the gas expansion coefficient, in particular, appearance of the factors 

in Eqs. (|65p - ([67]) is also quite revealing. It is the result of non-pcrturbative account of the influence exerted by 
the basic flow upon flame disturbances. Needless to say that such effects cannot be captured in principle by models 
based on weak-nonlinearity assumptions. 

Our investigation was based on the on-shell description of flames, developed in Refs. [1, H, 0, 0|j a-nd extended to 
the of anchored flames in Sec. Ill Bl This formulation allowed us to elucidate the role of the anchoring system and its 
influence on the flame structure, as well as to identify relevant boundary conditions for the flow variables. The simple 
and natural way this analysis was accomplished clearly demonstrates the power of this approach, not saying about 
that it permitted the analysis to be carried out at all. 

Another important technical aspect of our work is the locality issue discussed in Sees. Ill C[ IIII B 21 IIII C 21 As we 
have seen, the requirement of locality of the rod influence on the flame structure appears in the steady case analysis as 
the consistency condition . On the one hand, this condition expresses the fact that the piecewise constant gas flows 
of the basic V-pattern satisfy the main integro-differential equation p7p . and on the other hand, it serves for selecting 
inner solutions compatible with the given global flame structure. It is remarkable that the rod influence remains local 
also in the presence of flame disturbances. Namely, it was proved in Sec. IIII Cl that jumps in the functions u(x) and 
E(x) at a; = 0, which are potential sources of nonlocality, vanish in the high-velocity limit. 

The last important point to discuss is the practical conditions for applicability of the results obtained within the 
large-JJ limit. As is evident from the derivations of Sec. IIII C 11 in practical terms the condition V ^ 00 means 
that V should be large compared to [B — 1). At the same time, it is to be noted that validity of the asymptotic 
expansion of IK, obtained in Sec. IIIIB II requires only that U be large in comparison with unity. The latter condition 
is considerably weaker, taking into account that for real flames d is normally 5 to 8. This fact opens a way for 
investigation of moderate stream- velocities, which is the subject of the subsequent paper [l^. Another important 
issue is the influence of gravity. Recent experiments with open flames [l6l . \v\ demonstrate that the development of 
flame disturbances is strongly affected by the gravitational field. This effect can also be studied within our approach. 
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APPENDIX A: THE DARRIEUS-LANDAU RELATION 



In this appendix, we will demonstrate convenience of using two different imaginary units simultaneously for carrying 
out actual calculations. Namely, we will reproduce the classical result of linear stability analysis for planar flames, 
which will also serve as an important check of calculations that led us to Eq. ([35|) . 

In the case of freely propagating planar flames, one has s = 0, C/ = 1, ljo — so that Eq. ([35]) simplifies to 




I cot (^) + - V)] - 2i/'(a^)[ = , (Al) 



where H is the ordinary Hilbert operator, 

H exp{ikx) = ix{k) exp(ikx) , (A2) 

and we took into account the contribution due to {x — > —x)* by extending the range of ry-integration and doubling 
the last term. The linearized evolution equation takes the form 

uix) = iyf{x) . (A3) 

As usual, it is most convenient to look for a solution of these equations in a complex form. In doing so, however, 
one should be careful in respecting the original complex structure of Eq. (jAip . In order to preserve it, one can 
proceed in three different ways. The first is to extract the real and imaginary parts of Eq. (jAl[) . and then proceed to 
solving the system of equations in the usual way. This is the least convenient means, because it destroys the natural 
complex structure of Eq. (jAip . Another way followed in Ref . is to keep all intermediate relations involving the flow 
variables in an explicitly real form, like for instance in Eq. (jASp . The third method we choose here is to introduce a 
new imaginary unit, j, such that 

■2 -r •* ■ 

where the asterisk denotes the complex conjugation with respect to the initial imaginary unit, z, which had been used 
in the derivation of Eq. (jAlj) . 

I = — z, 

while the product (ij) is left unspecified. Thus, we write 

v = i'i+ji'2, u{x) — ue-'^'^ , u — ui+ju2, etc., 
where k is the wavenumber of perturbation, which according to the 2-periodicity condition takes on the values 

k = nm, m € Z. 

The physical solution is eventually found by extracting the real (or imaginary) part of the complex solution with 
respect to the unit j. 
One has 



+1 

I' 



rf^e^fe'Je-''^''/^ [zcot (^) + x{x - v)] = 77: \icot f ^"j [e^'^"*'^/^ - 5-^'=+"'/'^ 

L \&/ J jk — iv/0 L \0J 

-1 

}c)pi^jk~iu I9)x 
^ fr 
jk — iv jB ' 
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where the constant terms in the curly brackets cancel by virtue of the condition e^-''^ = 1. Using this in Eq. (|Aip yields 

2^' + - 1) (l + ^H) I '^'''^'^^"^ V = ■ 

V / 1 jkO — iv I 

Multiplying this equation by {jkO — iv), and extracting its real (with respect to i) part, we find 

2jkeu + 2ui~d' + (6* - 1) \viv{x) + jk9Hf'{x)^' = , (A4) 

while extraction of the imaginary part gives a similar equation, and comparison of the two leads to the relation 

w = Hu , 

which can be obtained also directly from (1 — iH)uj' = 0. Finally, writing u' ~ jku, /' — jkf, and expressing gas 
velocity via / with the help of Eq. (|A3p leads, after dividing by 0\k\f, to an algebraic equation 

^-tli,^ + 2iy\k\ - {0-l)k^ = O, 



from which the well-known Darrieus-Landau dispersion relation for the perturbation growth rate follows [id . 



i + e-l±i] \k\. 



APPENDIX B: EXTENSION OF EQ. dUD TO DISCONTINUOUS FUNCTIONS 

If the function a{x) in Eq. (|42|) does not satisfy conditions 

a(0+) = a(0-), a(+l) = a(-l), (Bl) 

its derivative is singular at x = 0, ±1, and the integration by parts used in the transition from Eq. (j4ip to Eq. (|42p is 
ambiguous. To correctly evaluate the integral, one has to turn back to the exact formula ([4]) in which all the functions 
involved are smooth, and apply it to a function A{x) satisfying (jBip . whose behavior near the rod or channel walls 
looks discontinuous from the outer point of view. More precisely, A{x) is supposed to vary rapidly for \x\ < R <^ I 
and near the walls, but normally at the intervals R < x < 1 — R and — 1 + i? < a; < —R, where it coincides with a{x). 
Thus, 

lim A(x) = a(x) . 

We also replace the function s|a;| describing the basic V-pattern by a smooth function F{x) such that 

lim F(a;) = slxl , a; € (-1, 1) . 

Neglecting the anchor dimensions means that the action of !K on a' is defined as 



^Ca') (x) ^ ]im{{^lA')] {x) . 



To find out how J-C acts on the derivative of A{x), we replace a by A in Eq. (|41]l . and integrate the right hand side by 
parts 

+1 

(^lA')ix) = lllpEl jd^A'{rj)cot[^{rj-x + t[F{rj)~F{x)])} 

-1 

+1 

= ~jd^[l + iF'{r,)]A{f^)cot[^{Tj-x + i[F{7j)-F{x)])}. (B2) 
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The boundary terms vanish here because the integral kernel is 2-periodic, and A{x) satisfies ^(—1) = A{+1), by 
the assumption. Since the functions A{x) and F'{x) have only finite jumps in the limit R —^ 0, the last integral in 
Eq. (|B2p is well-defined in this limit, representing a continuously difFerentiable function for all |a;| £ (0, 1). Thus, 



^hn { {^CA') Yx) = ~jd7^[l + isx{v)Hil) cot { I [77 



X + is{\ri\ 



Next, we go over to the large-slope limit. The right hand side of the last equation can be evaluated in this case in 
exactly the same way as we arrived to Eq. (|4T|) . Comparison with Eq. ([39|) shows that the role of the function 0(77) 
in this equation is now played by [1 -I- *sx('7)]a(?7), the only difference being that the large factor s comes from the 
integrand, rather than from the pre-integral factor in Eq. ((4]). Taking this into account, we readily find 



-za(-a;)(2|a;| - 1) 



d-q {a{ri)[sx{ri) - i] + a{-vi)[sx{-ri) - i]} xiv - 



-sx{x) {a{\x\) ~ a{-\x\)} + ix{x) {a{\x\) + a{~\x\)} 
-2ia{~x)x{x) + ia'{-x){2\x\ - 1) . 
Using the obvious identity x(a;){a(|a;|) + a(— — 2a(— a;)} = a(|a;|) — a(— we finally obtain 



{Via'^ (x) = {sx{x) - i) {a{-\x\) ~ a{\x\)\ + m'(-a;)(2|a;| - 1) 



which is exactly Eq. as was to be proved. Note that this result is independent of the particular choice of the 

functions A(x),F(x). 
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